A possible cooling effect in high temperature superconductors 
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We show that an adiabatic increase of the supercurrent 
along a superconductor with lines of nodes of the order pa- 
rameter on the Fermi surface can result in a cooling effect. 
The maximum cooling occurs if the supercurrent increases 
up to its critical value. The effect can also be observed in 
a mixed state of a bulk sample. An estimate of the energy 
dissipation shows that substantial cooling can be performed 
during a reasonable time even in the microkelvin regime. 



Although the mechanism of high temperature super- 
conductivity still remains unclear there is strong evidence 
that the order parameter in high temperature supercon- 
ductors has nodes on the Fermi surface |l| . Such a prop- 
erty is also peculiar to some superconductors with heavy 
fermions. In this paper we show that presence of nodes 
of the order parameter can lead to a cooling effect. The 
cooling is reached by adiabatically increasing the super- 
current around a superconducting ring or a cylinder. 

Currently the lowest temperature of a solid (T ~ ^I^K) 
is achieved by using a method of adiabatic nuclear de- 
magnetization [^U^. Spontaneous magnetic ordering of 
the nuclear magnetic moments represents the low tem- 
perature limit for nuclear refrigeration. The ordering 
temperature due to nuclear dipole-dipole interaction is 
typically a fraction of a microkelvin and, therefore, to 
achieve temperatures lower than O.lfiK a new refriger- 
ation technology is needed. In our cooling mechanism 
the conduction electrons are cooled during the adiabatic 
increase of the supercurrent instead of the nuclear spin 
system being cooled in the nuclear demagnetization pro- 
cess. Further development of the idea of cooling using the 
phenomena of superconductivity (superfluidity) could be 
promising to achieve the lowest temperature of solids. 

Cooling effect in clean superconductors 

Let us consider a superconducting ring (or a cylinder) 
made from a clean high temperature superconductor. For 
estimates one can assume the Fermi surface to be a cylin- 
der and the order parameter has a simple d-wave form 
A(p) = Aq cos{2(j>), where cj) is the polar angle in the ab 
crystalline plane. Taking into account the real form of 
the Fermi surface might change numerical coefficients, 
but these details are not necessary for estimating the 
magnitude of the cooling effect. Let the c -axis of the 
superconductor be parallel to the symmetry axis of the 
ring and the temperature of the system is T ^ Tc, where 
Tc is the superconducting transition temperature. 



At low temperatures the main contribution to the elec- 
tronic specific heat of the superconductor arises from 
quasiparticles that are activated above the gap in nar- 
row vicinities of the nodes of the order parameter. The 
energy of these quasiparticles is given by 



E{p) - + |A(p)i^ 



and the electronic specific heat is 



C(T) = 36C(3)Ar^ 
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(1) 
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where Np is the density of states at the Fermi surface 
for normal metal and C(3) — 1.202 is the Riemann zeta 
function (throughout this paper, apart from formula ([2l|), 
we use units in which = 1)- Therefore, the entropy of 
the system has the form: 



SiT) = f 
Jo 



C{T) 
T 



dT 



18C(3)7Vf^. 
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Let us now consider the same superconducting ring 
with a supercurrent uniformly distributed across the 
cross section of the ring. This situation is possible if 
the thickness of the ring d is less than penetration depth 
of the magnetic field A: c? < A. In the presence of such 
superflow with a superfluid velocity Vg, the quasiparticles 
can be considered quasiclassically with the energy 



E{p, Vg) = VC(P)' + |A(P)|2 + . 



(4) 



where Pp^_ is the component of the Fermi momentum 
in the direction perpendicular to the c-axis. Due to the 
existence of the supercurrent, quasiparticle states with 
negative energy appear for momentum directions within 
the vicinity of nodes of the order parameter. These states 
are occupied even at zero temperature. The presence of 
these states gives a term, linear in T, in the low tem- 
erature specific heat, which is dominant if pf±.Vs 2> T 
After a space averaging over the ring this term can 
be written in the form C{T) = 'i-n-NpTpF^Vs/SAo. As a 
result, the entropy in the presence of the supercurrent is: 



S[T) = -ttNp— — T. 



(5) 



Let us assume that at an initial moment of time there 
is no superflow along the ring, and the ring is thermally 
isolated. A magnetic field is then turned on adiabatically 
so that the magnetic fiux through the ring $ slowly in- 
creases in time. The induced electric field E ~ <^/lc {I is 
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the ring's length) generates current inside the ring which 
consists, in general, of superconducting and dissipative 
parts. The dissipative part of the current Idiss is propor- 
tional to E, and the energy dissipation can be estimated 
as 

/ I^ssEldt^ / $2^t ^ -1. (6) 
Jt=o Jt=o 

The dissipation can be made negligibly small if the time 
tf to build up the magnetic field is big enough. In this 
process all the work which is produced by the induced 
electric field goes to increase the supercurrent around the 
ring; and the part of this work that dissipates into heat 
and changes the entropy can be made negligibly small. 
The necessary cooling time tf is estimated at the last 
section of this article. 

The equation of entropy conservation in the adiabatic 
process of increasing the supercurrent has the form: 

ISCmF^ = ^nN.P^Tf, (7) 

Ao O Ao 

where Tq and Tf are the initial and the final temperatures 
of the ring. As a result, the final temperature Tf after 
the adiabatic increase of the supercurrent is 

This formula is applicable iipF±Vs ^ Tf. The maximum 
cooling occurs when the supercurrent (and the supcrfluid 
velocity) is equal to its critical value, that is pfi_Vs ^ 
2Ao ^ 4.28Tc and 

T/„.i„ « ^ (9) 

So, to achieve substantial cooling, it is important to use 
a superconductor with high Tc. If, for example, Tq — 
IK and Tc = lOOiiT, after the adiabatic increase of the 
supercurrent, the ring cools down to Tf — Q.Q\K. 

One should mention that we base our estimates on the 
assumption of uniform temperature distribution during 
the cooling process. However, the weight of quasiparticle 
states shifted from above the Fermi level to below de- 
pends on the angle the supercurrent makes with the nodal 
directions, which changes around the ring. This would 
provide a larger cooling effect at four points on the ring 
(at which Vs constitutes 45° angle with respect to a, h 
axes) compared to elsewhere and in the ideal case of zero 
thermoconductivity the temperature distribution over 
the ring would be nonuniform with T„iax/T„im — V^. 
However, finite thermoconductivity equalizes the temper- 
ature around the ring and nonuniform cooling would be 
very difficult to observe experimentally. Estimates show 
that nonuniform cooling can be detected if the cooling 
time is less than about R'^/vpTN, where R is the radius 
of the ring, vf is the Fermi velocity and tn is the relax- 
ation time in the normal state. For reasonable parame- 
ters we obtain that the necessary cooling time should be 



of the order of nanosecond. Real cooling time is much 
larger than nanosecond which guarantees approximately 
uniform temperature distribution during the cooling pro- 
cess. As a result the entropy dissipation due to the ther- 
mal current is negligible and does not modify our formu- 
las. 

Also we do not consider here the effect of fluxoid quan- 
tization which is negligible when the flux through the 
ring is much greater than the flux quantum. This gives 
a restriction on the superfluid velocity ^ vfci/R, 
where a is the interatomic spacing. From the other 
hand, the critical superfluid velocity is of the order of 
2Ao/pf± ~ vfcl/S,, where ^ is the coherence length. So, 
we can omit the effect of fluxoid quantization in our prob- 
lem if R > ^. 

One should also note that value Vg ^ 2Aq/pf± de- 
termines depairing current density jdepair- Because of 
the grain structure of high temperature superconductors 
the critical current density jc is usually much less than 
the depairing value. A misalignment of adjacent grains 
suppresses the order parameter at the grain boundary 
and results in Josephson tunneling behavior and a near- 
exponential decrease in jc with misorientation angle. If 
the current through the superconductor is increased to jc, 
rather than jdepair , the estimate (^) for the cooling tem- 
perature should be multiplied by the factor jdepair/ jc- To 
reach jc close to the depairing value one should use ma- 
terials with well- aligned grains, which can be achieved in 
thin films For example, YBa2Cu30r-s films were ob- 
tained with the critical current density jc ~ 40MA/cm^, 
so that jdepair/ jc ~ 5 |^ . It is worth to note that if the 
cooling is achieved by applying external magnetic field, 
rather than creating a supercurrent along the film, the 
grain structure would not be an obstacle. In such a way, 
the magnetic field induces screening supercurrent in each 
grain. The intragrain supercurrent is not limited by the 
small value of the Josephson critical current in weak links 
between the grains and can be close to jdepair ■ 

It is interesting to estimate the value of the cooling 
effect for a bulk sample in a mixed state. In a mag- 
netic field H {Hci,Hc2T^/T^ < i/ < i?c2) the en- 
tropy of a clean superconductor with line(s) of nodes is 
S cx TyjH/Hc2 ||]. So, an adiabatic increase of mag- 
netic field from Hq to Hf cools down the sample to the 
temperature 




One should note that in conventional s -wave supercon- 
ductors the electronic contribution to the entropy is ex- 
ponentially small (due to finite gap in the excitation spec- 
trum) and thus not dominant at low temperatures. How- 
ever, in a mixed state of s -wave superconductors normal 
regions in the vortex cores give the linear temperature 
contribution to the electron entropy S ^ NfTH/Hc2, 
which can exceed the phonon contribution at low enough 
temperature T (typically less then IK) 
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Td 
Tf 



1/2 



H 

Hc2 



1/2 
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where To and Tp are Debye and Fermi temperatures 
accordingly. In this temperature region an adiabatic in- 
crease of apphed magnetic field H {Hd < H < Hc2) 
results in the cooling effect Tf = ToHo/Hf. 

Effect of impurities 

At very low temperatures (when T < 7) the effect of 
impurity-induced bound states becomes significant. Let 
us consider the cooling effect at temperatures when the 
regime T < 7 <g; is achieved. The energy scale 7 
is the bandwidth of quasiparticle states bound to im- 
purities P,p^. For an order parameter with a line of 
nodes, the bandwidth 7 and density of the quasiparti- 
cle states at zero energy A^(0) are finite for any nonzero 
concentration of impurities. The quantities 7 and A^(0) 
depend on the impurity concentration and the scat- 
tering phase shift 60 . For example, for scattering in the 
unitary limit (^o — 7i'/2) and the d-wave order parame- 
ter, 7 ~ V^r«An/2 , N{0) = 8-fNF\og{AAo/^) /nAo, 
where = Ui/nNp. In the opposite limit (Born ap- 
proximation with So — > 0) the bandwidth 7 is exponen- 
tially small 7 ^ 4Ao exp(— 7rAo/2F„^g) and, therefore, 
the regime T < 7 can be very difficult to achieve experi- 
mentally. For most high temperature and heavy fermion 
superconductors the impurity scattering is strong (or 
an intermediate strength) rather than in the Born limit 

0-0- 

For T < 7 and in the absence of superfiow the elec- 
tronic specific heat is given by C(T) = tt^N{0)T/3. So, 
instead of (||), the entropy of the system is 



S{T) = —N{0)T. 



(12) 



After the adiabatic process of creating the supercurrent 
around the ring with pfi_Vs 3> 7, the entropy of the 
system is given approximately by the same expression 
(g) as for the clean superconductor. As a result, the 
equation of the entropy conservation gives the following 
formula for the final temperature: 



7r7V(0) Ao ^ 

If — Jo 

4 Nf Pf±Vs 
For maximum cooling (jjf_lVs ^ 2Ao) 

7riV(0) 



f min 



Nf 



-Tn 



For scattering in the unitary limit we get: 



_ 7 /4Ao 
T,^i„«.-log(^ — 



To 



T, 



■To 



(13) 



(14) 



(15) 



In the case of a bulk superconductor in a mixed state 
the presence of impurities results in appearance of a 



crossover field H* when square-root magnetic field de- 
pendence of the specific heat changes to a dependence of 
the form HlnH |l5|. U H > H* the influence of impu- 
rities is negligible and the cooling effect is given by Eq. 
(|l^). However, if iJo < H* < Hf then Tf « Toy/H*/Hf. 

Estimation of the energy dissipation 

If there is a time-dependent current in a superconduc- 
tor, the normal electrons give a finite amount of dissipa- 
tion because the supercurrent is not a zero-impedance 
shunt at nonzero frequency oj (in a two-fluid model 
Ima — UsE^ jmuS). This dissipation increases the en- 
tropy of the system. Given an imposed ac current den- 
sity J = jo sin (wi), the power dissipated per unit volume 
at low frequencies is |l6||: 



dW, 



diss 



dt 



Rea 9 




2\4 



(16) 



where a is the conductivity of the superconductor, A is 
the penetration depth. To estimate the change of the 
entropy AS* in the process of increasing the supercur- 
rent from j — to j — jo one can integrate Eq. ( [l^ ) 
in the limits < t < tf = 7r/2w. Generally speaking, 
the conductivity a can depend on the magnitude of the 
supercurrent in the superconductor. However for estima- 
tion one can take a to be equal to its maximum value (at 
a maximum value of the current). As a result, we obtain 
the following estimate for the change of the entropy per 
unit volume: 

As^r^'^dt^ 



T dt 



leTT^A" 



sin (ujt) Reafdt « ^XftfC* 

(17) 

The energy dissipation gives rise to a heating of the su- 
perconductor and changes the final temperature Tf. To 
estimate the magnitude of this effect one should add 
AS* to the left side of Eq. (^. It is possible to ob- 
tain substantial cooling even if AS* is larger than the ini- 
tial entropy. However, to make the cooling effective, one 
should increase the supercurrent slowly enough to ensure 
that the dissipation does not significantly change the en- 
tropy: AS < 4:TrNFTfPF±Vs/3Ao. Taking into account 
jo ~ ensVs — mc^Vs/ AneX^ , we see that dissipation can 
be omitted if the time of increasing the supercurrent from 
zero to its critical value satisfies the following condition 



Rea 



^0 



(18) 



The electrical conductivity for a superconductor with an 
order parameter that vanishes along a line of nodes has a 
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universal limiting value which is independent of the con- 
centration of impurities and the scattering phase shift 
pTf . For the two dimensional d-wave order parameter 
and an isotropic 2D Fermi surface the universal value 
has the form Rea (w 0, T < 7) ~ c^NfvIta, 

where ta — H/ttAq is a universal transport time. How- 
ever, the superflow modifies this result and necessarily 
increases the conductivity in the region interesting for 
cooling (i.e. when pf±Vs ^ 7)- 

To estimate the electrical conductivity in the presence 
of the superflow, one can use the same semiclassical ap- 
proach as for the case of thermal conductivity At 
uj —^ and 7 <^ Pfj^Vs the final result is given by (cf. 

El): 



Reai 



2mT 



doj- 



cosh 



\2T J 



2tt 



Plpj 



\uj - Pfi_Vs 



{Co ~Pfi_Vs 



(19) 



where r(a), u^) is the relaxation time for quasiparti- 
cle scattering on impurities in the superconducting 
state, n is the electron density. For the d-wave or- 
der parameter we obtain for T <C Pf±Vs ^ Aq, 
Reau = ne'^Tpf /tn, in the Born approximation; Recu = 
ne'^Tpf {1 + \ cos pjpj_v'^ / Am Aq, in the unitary limit, 
where tn is the relaxation time in the normal state, 
is the angle between Vg and the i-axis, i = a, b. 
For the supercurrent close to its critical value we get 
the estimate Rea ^ ne^TN/ra = e'^NFVpTN which is 
tn/ta ^ Aq/7^ ^ 1 times larger than the universal 
value. If we substitute this estimate into (p^), we obtain 
the following limitation for the necessary cooling time 



tf > tnT^JT]. 

In the unitary limit r^r = 7rfi.Ao/4A:B7^ 
the limitation can be rewritten as 



> 



t2 hn 



r 



(20) 

hTc/ks'j'^ and 



(21) 



If, for example, ~ lOOK, 7 = O.OlTc, Tq = O.OIK, 
Tf ~ 10~'^K the estimate for the necessary cooling time 
is tf > 12min. In principle, the cooling time can be rea- 
sonable even in a submicrokelvin region. In this region 
unconventional superconductors with low Tc like heavy 
fermion systems can be used. One of the possible candi- 
dates is UPt^, which has an order parameter with a line 
of nodes and Tc ~ 0.5K. If, for example, Tc = 0.5K, 
7 = O.lTc, To = 10~^K, Tf - 10-"^ K then from Eq. 
( pl| ) we obtain t/ > lO/i. One should note, however, that 
applicability of the method is limited by ultra low tem- 
peratures. At microkelvin temperatures the entropy of 
the quasiparticles near nodes is small as compared to the 
nuclear entropies and an estimate of the Korringa con- 
stant for UPt^ from the value of pure Pt (K = 3QmKs) 



shows that in 10 hours the electrons are not thermally 
disconnected from the nuclear system. 

To make effective refrigerating storage based on the 
method discussed here, it is not sufRcient to use a single 
ring or a cylinder because the allowed width of the ring 
is confined by the penetration depth A, which is usually 
a fraction of a micron. Instead of a single ring, one can 
make a coil from a superconducting wire or a tape so 
that the width of the tape is of the order of A, but the 
total width of the coil can be made much larger. The 
superconducting current in the tape can be created by 
applying an external magnetic field to the coil or by ap- 
plying voltage to the tape ends. Another possibility to 
obtain a cooling effect is to use a bulk sample in a mixed 
state. 

In conclusion, we have shown that an adiabatic in- 
crease of the supercurrent through a superconductor with 
an order parameter with line(s) of nodes on the Fermi 
surface results in a cooling effect. For maximum cooling, 
the initial temperature Tq decreases to Tq /Tc for clean 
superconductors and to Tq^/Tc for Tq < 7. The effect 
can also occur in a mixed state of a bulk sample. To 
observe the cooling effect it is important to be in the low 
temperature region in which conduction electrons give 
the main contribution to the system's entropy. 
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